Determinants of SARS-CoV-2 receptor gene expression in upper and lower airways (preprint)

The recent outbreak of the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), which causes coronavirus disease 2019 (COVID-19), has led to a worldwide pandemic. One week after initial symptoms develop, a subset of patients progresses to severe disease, with high mortality and limited treatment options. To design novel interventions aimed at preventing spread of the virus and reducing progression to severe disease, detailed knowledge of the cell types and regulating factors driving cellular entry is urgently needed. Here we assess the expression patterns in genes required for COVID-19 entry into cells and replication, and their regulation by genetic, epigenetic and environmental factors, throughout the respiratory tract using samples collected from the upper (nasal) and lower airways (bronchi). Matched samples from the upper and lower airways show a clear increased expression of these genes in the nose compared to the bronchi and parenchyma. Cellular deconvolution indicates a clear association of these genes with the proportion of secretory epithelial cells. Smoking status was found to increase the majority of COVID-19 related genes including ACE2 and TMPRSS2 but only in the lower airways, which was associated with a significant increase in the predicted proportion of goblet cells in bronchial samples of current smokers. Both acute and second hand smoke were found to increase ACE2 expression in the bronchus. Inhaled corticosteroids decrease ACE2 expression in the lower airways. No significant effect of genetics on ACE2 expression was observed, but a strong association of DNA- methylation with ACE2 and TMPRSS2- mRNA expression was identified in the bronchus.

Prime-boost protein subunit vaccines against SARS-CoV-2 are highly immunogenic in mice and macaques (preprint)

SARS-CoV-2 vaccines are advancing into human clinical trials, with emphasis on eliciting high titres of neutralising antibodies against the viral spike (S). However, the merits of broadly targeting S versus focusing antibody onto the smaller receptor binding domain (RBD) are unclear. Here we assessed prototypic S and RBD subunit vaccines in homologous or heterologous prime-boost regimens in mice and non-human primates. We find S is highly immunogenic in mice, while the comparatively poor immunogenicity of RBD was associated with limiting germinal centre and T follicular helper cell activity. Boosting S-primed mice with either S or RBD significantly augmented neutralising titres, with RBD-focussing driving moderate improvement in serum neutralisation. In contrast, both S and RBD vaccines were comparably immunogenic in macaques, eliciting serological neutralising activity that generally exceed levels in convalescent humans. These studies confirm recombinant S proteins as promising vaccine candidates and highlight multiple pathways to achieving potent serological neutralisation.

An intestinal cell type in zebrafish is the nexus for the SARS-CoV-2 receptor and the Renin-Angiotensin-Aldosterone System that contributes to COVID-19 comorbidities (preprint)

People with underlying conditions, including hypertension, obesity, and diabetes, are especially susceptible to negative outcomes after infection with the coronavirus SARS-CoV-2. These COVID-19 comorbidities are exacerbated by the Renin-Angiotensin-Aldosterone System (RAAS), which normally protects from rapidly dropping blood pressure or dehydration via the peptide Angiotensin II (Ang II) produced by the enzyme Ace. The Ace paralog Ace2 degrades Ang II, thus counteracting its chronic effects. Ace2 is also the SARS-CoV-2 receptor. Ace, the coronavirus, and COVID-19 comorbidities all regulate Ace2, but we dont yet understand how. To exploit zebrafish (Danio rerio) as a disease model to understand mechanisms regulating the RAAS and its relationship to COVID-19 comorbidities, we must first identify zebrafish orthologs and co-orthologs of human RAAS genes, and second, understand where and when these genes are expressed in specific cells in zebrafish development. To achieve these goals, we conducted genomic analyses and investigated single cell transcriptomes. Results showed that most human RAAS genes have an ortholog in zebrafish and some have two or more co-orthologs. Results further identified a specific intestinal cell type in zebrafish larvae as the site of expression for key RAAS components, including Ace, Ace2, the coronavirus co-receptor Slc6a19, and the Angiotensin-related peptide cleaving enzymes Anpep and Enpep. Results also identified specific vascular cell subtypes as expressing Ang II receptors, apelin, and apelin receptor genes. These results identify specific genes and cell types to exploit zebrafish as a disease model for understanding the mechanisms leading to COVID-19 comorbidities. SUMMARY STATEMENTGenomic analyses identify zebrafish orthologs of the Renin-Angiotensin-Aldosterone System that contribute to COVID-19 comorbidities and single-cell transcriptomics show that they act in a specialized intestinal cell type.

Structural Variants in SARS-CoV-2 Occur at Template-Switching Hotspots (preprint)

The evolutionary dynamics of SARS-CoV-2 have been carefully monitored since the COVID-19 pandemic began in December 2019, however, analysis has focused primarily on single nucleotide polymorphisms and largely ignored the role of structural variants (SVs) as well as recombination in SARS-CoV-2 evolution. Using sequences from the GISAID database, we catalogue over 100 insertions and deletions in the SARS-CoV-2 consensus sequences. We hypothesize that these indels are artifacts of imperfect homologous recombination between SARS-CoV-2 replicates, and provide four independent pieces of evidence. (1) The SVs from the GISAID consensus sequences are clustered at specific regions of the genome. (2) These regions are also enriched for 5 and 3 breakpoints in the transcription regulatory site (TRS) independent transcriptome, presumably sites of RNA-dependent RNA polymerase (RdRp) template-switching. (3) Within raw reads, these structural variant hotspots have cases of both high intra-host heterogeneity and intra-host homogeneity, suggesting that these structural variants are both consequences of de novo recombination events within a host and artifacts of previous recombination. (4) Within the RNA secondary structure, the indels occur in "arms" of the predicted folded RNA, suggesting that secondary structure may be a mechanism for TRS-independent template-switching in SARS-CoV-2 or other coronaviruses. These insights into the relationship between structural variation and recombination in SARS-CoV-2 can improve our reconstructions of the SARS-CoV-2 evolutionary history as well as our understanding of the process of RdRp template-switching in RNA viruses.

Modeling SARS-CoV-2 infection in vitro with a human intestine-on-chip device (preprint)

Coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus (SARS-CoV-2) has given rise to a global pandemic. The gastrointestinal symptoms of some COVID-19 patients are underestimated. There is an urgent need to develop physiologically relevant model that can accurately reflect human response to viral infection. Here, we report the creation of a biomimetic human intestine infection model on a chip system that allows to recapitulate the intestinal injury and immune response induced by SARS-CoV-2, for the first time. The microengineered intestine-on-chip device contains human intestinal epithelium (co-cultured human intestinal epithelial Caco-2 cells and mucin secreting HT-29 cells) lined in upper channel and vascular endothelium (human umbilical vein endothelial cells, HUVECs) in a parallel lower channel under fluidic flow condition, sandwiched by a porous PDMS membrane coated with extracellular matrix (ECM). At day 3 post-infection of SARS-CoV-2, the intestine epithelium showed high susceptibility to viral infection and obvious morphological changes with destruction of intestinal villus, dispersed distribution of mucus secreting cells and reduced expression of tight junction (E-cadherin), indicating the destruction of mucous layer and the integrity of intestinal barrier caused by virus. Moreover, the endothelium exhibited abnormal cell morphology with disrupted expression of adherent junction protein (VE-cadherin). Transcriptional analysis revealed the abnormal RNA and protein metabolism, as well as activated immune responses in both epithelial and endothelial cells after viral infection (e.g., up-regulated cytokine genes, TNF signaling and NF-kappa B signaling-related genes). This bioengineered in vitro model system can mirror the human relevant pathophysiology and response to viral infection at the organ level, which is not possible in existing in vitro culture systems. It may provide a promising tool to accelerate our understanding of COVID-19 and devising novel therapies.

Spike mutation D614G alters SARS-CoV-2 fitness and neutralization susceptibility (preprint)

A spike protein mutation D614G became dominant in SARS-CoV-2 during the COVID-19 pandemic. However, the mutational impact on viral spread and vaccine efficacy remains to be defined. Here we engineer the D614G mutation in the SARS-CoV-2 USA-WA1/2020 strain and characterize its effect on viral replication, pathogenesis, and antibody neutralization. The D614G mutation significantly enhances SARS-CoV-2 replication on human lung epithelial cells and primary human airway tissues, through an improved infectivity of virions with the spike receptor-binding domain in an "up" conformation for binding to ACE2 receptor. Hamsters infected with D614 or G614 variants developed similar levels of weight loss. However, the G614 virus produced higher infectious titers in the nasal washes and trachea, but not lungs, than the D614 virus. The hamster results confirm clinical evidence that the D614G mutation enhances viral loads in the upper respiratory tract of COVID-19 patients and may increases transmission. For antibody neutralization, sera from D614 virus-infected hamsters consistently exhibit higher neutralization titers against G614 virus than those against D614 virus, indicating that (i) the mutation may not reduce the ability of vaccines in clinical trials to protect against COVID-19 and (ii) therapeutic antibodies should be tested against the circulating G614 virus before clinical development. ImportanceUnderstanding the evolution of SARS-CoV-2 during the COVID-19 pandemic is essential for disease control and prevention. A spike protein mutation D614G emerged and became dominant soon after the pandemic started. By engineering the D614G mutation into an authentic wild-type SARS-CoV-2 strain, we demonstrate the importance of this mutation to (i) enhanced viral replication on human lung epithelial cells and primary human airway tissues, (ii) improved viral fitness in the upper airway of infected hamsters, and (iii) increased susceptibility to neutralization. Together with clinical findings, our work underscores the importance of this mutation in viral spread, vaccine efficacy, and antibody therapy.